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and 8° to the north of the Bar. Observed metaUicity distributions show a well 
defined peak, with a tail toward low metallicities. The mean metaUicity remains 
constant until 6° ([Fe/H]~-0.5 dex), while for the outermost field, at 8°, the 
mean metaUicity is substantially lower than in the rest of the disk ([Fe/H]~-0.8 
dex). The combination of spectroscopy with deep CCD photometry has allowed 
us to break the RGB age-metaUicity degeneracy and compute the ages for the 
objects observed spectroscopicaUy. The obtained age-metaUicity relationships 
for our four fields are statistically indistinguishable. We conclude that the lower 
mean metaUicity in the outermost field is a consequence of it having a lower frac- 
tion of intermediate-age stars, which are more metal-rich than the older stars. 
The disk age-metaUicity relationship is similar to that for clusters. However, 
the lack of objects with ages between 3 and 10 Gyr is not observed in the field 
population. Finally, we used data from the literature to derive consistently the 
age-metaUicity relationship of the bar. Simple chemical evolution models have 
been used to reproduce the observed age-metaUicity relationships with the pur- 
pose of investigating which mechanism has participated in the evolution of the 
disk and bar. We find that while the disk age-metaUicity relationship is well 
reproduced by close-box models or models with a small degree of outfiow, that 
of the bar is only reproduced by models with combination of infall and outfiow. 

Subject headings: galaxies: evolution — galaxies: stellar content — Magellanic 
clouds 



1. Introduction 



Despite decades of work, there are still signifi cant gaps in our know ledge of the LMC's 
star formation and chemical enrichment histories (jOlszewski et al.lll996l ). This is motivated 
in part by the vastness of its stellar populations and to our limitations in observing siz- 
able samples of st ars. The age distribution of star clusters is relatively well known (e.g. 
Geisler et al.l 119971 ): there is an age interval between 10 and 3 Gyr with almost no clusters . 
This so-called age-gap may also correspond to an abundance gap (jOlszewski et al.lll99ll ). 
since the old clusters are metal-poor while the young ones are relatively metal-rich. The star 
formation history (SFH) of field stars is much less precisely known. Studies using HST data 
in small fields, suggest that the SFH of the LMC disk has been more or less continuous, with 



some increase in the star formation rate (SFR] in the last few Gyr (e.g. ISmecker-Hane et al. 
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Holtzman et al. 


1999) 



19991 ). This contrasts with previous results (based 



on much shallower data), which found a relatively young age (a few Gyr) for the dominant 
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LMC population flHardv et alJll984l : ISertelli et al.lll992l : IWesterlund. Linde. fc Lvngalll995 



Vallenari et al.lll996l ). The SFH of the LMC can now be studied using sufficiently deep 



ground-based data (e.g. reaching the oldest main- sequence turn-off wit h good photometric 
preci sion) in large areas and in different positions of the galaxy (e.g. 
20071 . hereafter Paper I). 



Gallart et al. 2004 



Less is known about the LMC chemical enrichment history. The ag e-metallicity re 



lation (AMR) norrnally used for the LMC is defined by star clusters (e.g. [Olszewski et al. 
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ters obtain similar results. The mean metallicity jumps from [Fe/H] ~ —1.5 for the old- 
est clusters to [Fe/H] ~ —0.5 for the y oungest ones. There were no studies on the field 



population AMR until the last decade. iDopita et al.l (Il997l ). from a study of a-elements 
in planetary nebulae, obtained a result qualitative ly similar t o tha t found in the clus- 
ters although on ly ten objects were used. Later, iBica et al.l (119981) . using Was hington 
photometry, and ICole. Smecker-Hane. fc Gallagher I (120001 ) and iDirsch et al.l (120001 ) , using 
Stron grem photometry, obtained the me t allicit y of RGB stars in different positions of the 
LMC. ICole. Smecker-Hane. &: Gallagher I (l2000l ) also observed 20 stars with the infrared Call 
triplet (CaT). All of them found t hat the age-gap ob served in the clusters is not found in 
the field population. More recently, ICole et al.l (120051 ) have obtained stellar metallicities for 
almost 400 stars in the bar of the LMC, also using CaT lines. The metallicity for each star 
has been combined with its position in the color-magnitude diagram (CMD) to estimate 
its age. The obtained AMR shows a similar behavior to that of the clusters for the oldest 
population. However, while for the cluster s the metallicity has i ncreased over the last 2 Gyr, 
this has not happened in the bar. Finally iGratton et al.l (120041 ) measured the metallicity of 
about 100 bar RR Lyrae, obtaining an average metallicity of [Fe/H] = —1.48. 

Another point that it is necessary to investigate is the p resence, or not, o f an abun- 
dance gradient i n the LMC. From observations of clusters, lOlszewski et al.l ( 1l99ll ) and 
Santos Jr. et al.l (119991 ) found no evidence for the presence of a radial metallicity gradient in 
the LMC. T he only evidence for a radial metallic i ty gra dient in the LMC cluster system was 
reporte d bvlKontizas. Kontizas. fc Michalitsiano J1993I ) based on six outer LMC clusters (> 
8 kpc). iHill. Andrievsky. &: Spitd (119951 ) were the ffist to report evidences of a gradient in 
the field population from high-resolution spectroscopy, using a sample of nine stars in the 
bar and in the disk. They found that the bar is on average 0.3 dex more metal-rich than 
the disk population, but t he di sk stars studied are located within a radius of 2°. Subsequent 
work by lCioni fc Habind (120031 ) . detected that the C/M ratio between number of asymptotic 
giant branch stars of spectral types C and M increased when moving away from the bar and 
within a radius of 6?7. As the C/M ratio is anticorrelated with metallicity, this increment 
implies a decrease in metallicity. Finally, an outward radial gradient of decreasing metallicity 
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was also found by lAlvea (120041 ) from infrared CMD using the 2MASS survey. 



We have obtained deep photometry with the Mosaic II CCD Imager on the CTIO 4m 
telescope in four disk fields at different distances from the center of the LMC (RA=5''23'"34!5, 
5=-69°45'22") with a quality similar to that obtained by HST in more crowded areas (see 
Paper I). The position of these fields, together with a description of their CMDs, are presented 
in Paper I, and detailed SFHs will be published in forthcoming papers. In the present 
investigation we focus on obtaining stellar metallicities for a significant number of individual 
RGB stars in these four fields using spectra obtained with the HYDRA spectrograph at the 
CTIO 4m telescope. In Section [2] we present our target selection. The observations and 
data reduction are presented in Section [31 The radial velocities of the stars in our sample 
are obtained in Section HI In Section O we discuss the calculation of the CaT equivalent 
widths and the determination of metallicities. Section [HI presents the method used to derive 
ages for each star by combining the information on their metallicity and position on the 
CMD. The analysis of the data is presented in Section [71 where the possible presence of 
a kinematically hot halo is discussed and the AMR and the SFH in each of our fields are 
described and compared. In the final section, the AMRs are compared to theoretical models 
and conclusions are drawn about the chemical evolution of the LMC. 



2. Target Selection 

Our starting point are deep CMDs of four 36' x 36' fields located ~3°, 5°, 6° and 8° 
from the LMC center, which are presented and discussed in Paper I. All fields are located 
north of the LMC bar, and their CMDs reach the oldest main-sequence turn-offs. A clear 
gradient exists in the amount of intermediate-age and young populations from the inner to 
the outermost fields, in the sense that the young population is more prominent in the inner 
parts (see Paper I for a more detailed discussion). 

In each field we have selected the stars to be observed with HYDRA in two windows in 
the CMD, which are plotted in Figure [H The bluest color limit has been selected to avoid 
upper red-clump and red supergiant stars and the lower magnitude limit to avoid including 
red-clump stars in the spectroscopic sample. The blue edge of the window is sufficiently 
blue to ensure that no metal poor stars would be excluded (the position of the blue ridge 
corresponds approximately with the position of ~2 Gyr old stars with [Fe/H]=-3). Similarly, 
the opposite edge is sufficiently red so that even old stars with solar metalicities, or even 
higher, fit within the box. These selected stars have been ordered from the brightest to the 
faintest ones, with no color restriction. This list has been used as input for the configuration 
task of the instrument, which tries to optimize the number of allocated fibers. The first stars 
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in the list, the brightest ones, have priority over the others. 



3. Observations and Data Reduction 

The observations were carried out in two different runs, in December 2002 and in January 
2005, at the CTIO 4 m telescope with the HYDRA multifiber spectrograph. We used the 
KPGLD grating, which provides a central wavelength of 8500 A, and a OG590 order-blocking 
filter. The physical pixels were binned 2 x 1 in the spectral direction yielding a dispersion 
of 0.9 A/pix. HYDRA has a field of view of about 40 arcmin, which on account of the high 
density of stars in our fields, allowed us to observe more than 100 stars in each configuration, 
with the 140 available fibers. Because of technical problems and bad weather, on our first 
run we were able to observe only two configurations located at 5° and 8° respectively. In 
the more successful second run, one fiber configuration in each field was observed except for 
the closest one, where we observed two configurations. Stars in two globular clusters, NGC 
4590 and NGC 3201, were also observed. Those in the first were used as radial velocity 
standards, but also to compare the equivalent widths determi ned through the pres ent set-up 



with those setups used for the cluster observations reported in lCarrera et al.l (120071 . hereafter 
Paper II). Since the differences were negligible (see Figure 1 of paper II), the second cluster 
was included in the calibration of the CaT as metallicity indicator. We concluded that 
the LMC equivalent widths obtained with the present observational setup could be directly 
transformed into [Fe/H] using the calibration in paper II. Equivalent widths, magnitudes 
and radial velocities for the 702 stars observed in the LMC are listed in Table [H 

The data reduction was done following the procedure described in Knut Hydra notej^. 
First, c osmic rays were rem oved from the images using the IRAF@ Laplacian edge-detection 
routine ( van Dokkunlbool ). All images were then bias- and overscan-subtracted, and trimmed. 



Flat exposures were taken during the day with a diffusing screen installed in the spectro- 
graph to obtain the so-called "milk-flat", which is used to correct the effect of the different 
sensitivities among pixels. It was applied to all images using ccdproc. The spectra were 
extracted, flat-field corrected (see bellow), and calibrated in wavelength using dohydra, a 
program developed specifically to reduce the data obtained with this spectrograph. The 
lamp flats acquired at the beginning of each configuration were used to define and trace each 



^http://www. ctio.noao.edu/spectrographs/hydra/hydra-knutnotcs. html 

^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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aperture. The sky flats were used to eliminate the residuals after the flat-field correction with 
the milk-flat. Arcs, obtained before and after each configuration, were used for wavelength 
calibration. 

Dohydra can also perform sky subtraction. However, we noticed that after sky subtrac- 
tion, there remained significant sky line residues. We therefore developed our own procedure 
to remove the contribution of the sky lines from the object spectra. For a given configuration, 
we obtain an average high signal-to-noise ratio (S/N) sky spectrum from all fibers placed on 
the sky. Before subtracting this averaged sky from each star spectrum, we need to know the 
relation between the intensity of the sky in each fiber (which varies from fiber to fiber owing 
to the different fiber responses) and the average sky. This relation is used as a weight (which 
may depend on wavelength) that multiplies the average sky spectrum before subtracting it 
from each star. Our task minimizes the sky line residuals over the whole spectral region 
considered and allows very accurate removal of the sky emission lines. An example of raw, 
sky and final object is shown in Figure [2l 

We obtained between three and four exposures of 2700 s in each configuration. To 
minimize the contribution of residuals arising from cosmic rays and bad pixels, all spectra 
of the same object were combined. Finally, each object spectrum was normalized with the 
continuum task by fitting a polynomial, making sure to exclude all absorption lines in our 
wavelength range, such as the CaT lines themselves. 



Radial Velocities 



Our main purpose in obtaining radial velocities for our program stars is to reject pos- 
sible foreground objects. However, radial velocities can provide more useful information, as 
discussed in depth in Section flA] The radial velocities wer e calculated with the p ccor task 
in IRAF, which performs the cross-correlation described by lTonry fc David (Il979l ) between 
the spectra of the target and selected templates of known radial velocities. As templates, we 
used three st a rs in the globular cluster NGC 4590 (M68) with radial velocities obtained by 
Geisler et al.l ( 1l995l ). The final radial velocities for all the stars in our sample were obtained 
as the average of the velocities obtained with the three templates, weighted by the width of 
the corresponding correlation peaks. Radial velocity histograms for each field are plotted in 
Figure O A Gaussian function was fitted in order to obtain the central peak and the a of 
each distribution (values listed in Table [2]). The velocity distribution of the field located 3° 
from the center shows two peaks, w hich are marked i n Figure [31 Stars with radial velocity 
in the range 170 <Vr< 380 km s~^ (jZhao et al.ll2003l ) are considered LMC members. Only 
about 20 stars in each configuration of the instrument have been excluded based on their 
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radial velocity. These results are discussed in depth in Section mi 



5. CaT Equivalent Widths and Metallicity Determination 



The metallicity of the RGB stars is obtained following the procedure described in Pa- 
per II. In short, the equivalent width is the line area normalized to the continuum. The 
continuum is calculated from the linear fit to the mean value of each bandpass, defined to 
ob tain the contiriuum p osition. We have used the continuum and line bandpasses defined 
by ICenarro et al.l (120011 ) , which are listed in Table [3l The equivalent widths are calculated 
from profile fitting using a Gaussian plus a Lorentzian. This combination provides the best 
fit to line core and wings as discussed in Paper II. The CaT index is defined as the sum of 
the equivalent widths of the three CaT lines, denoted as SCa. The SCa calculated for each 
observed star and its uncertainty are given in Table [H together with the star magnitude and 
radial velocity. The reduced equivalent width, Wj, for each star has been calculated as the 
value of SCa at M/=0, using the slope obtained in Paper II for the calibration clusters in 
the Mi-ECa plane = -0.611 Amag^^). 

In Paper II we obtain ed relationships between the reduced equivalent widths, Wy and 
W^, and metallicity , onthe lzinn fc WestI (|l984l ). ICarretta fc Grattonl (|l997l . hereafter CG97) 
and iKraft &: Ivand (120031 ) scales. In this case, we will only use the relationships obtained 
on the CG97 metallicity scale, because it is the only one in which the metallicities of our 
open and globular calibration clusters have been obtained in a homogeneous way from high- 
resolution spectroscopy. 

In order to verify that the V and I relationships give similar results, we have applied both 
to obtain the metallicities of the stars observed in the field located at 5°. The differences are 
on average —0.01 dex, and never larger than 0.15 dex for stars with {V — I) < 2.5. For stars 
with colors larger than this, the differences are on average 0.2 dex and never get above 0.35 
dex. The calibration stars used in Paper II have {V — I) < 2.5, so the metallicities calculated 
for the reddest stars are necessarily extrapolations of the relationships. On average, the 
relationship based on M/ yields metallicities slightly lower than those obtained from My, 
but always within the uncertainties. In what follows we will use the relationship based on 
M/, because: i) the photometric accuracy in I is slightly better than in V, especially for 
the reddest stars on the RGB; ii) the RGB is better resolved in the I band; and iii) from a 



theoretical point of view (jPont et al.l 120041 ). it is likely that the relationship based on M/ is 
less sensitive to age. In short, then, the metallicity for each star is given by: 
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[Fe/H]cG97 = -2.95 + 0.38SCa + 0.23Mj (1) 

Binarity could affect the derived metallicity only in the unlikely case of the two stars in 
the system being RGB stars of a very similar mass. In this case the system would be ~0.7 
mag brighter than a single star. As the binarity might not affect the equivalent width of the 
CaT lines, the derived metallicity would be lower by ~0.4 dex than the actual metallicity. 

Figure m shows the position of LMC radial velocity member stars in the Mj-SCa plane 
for our four fields, from the innermost one (top left) to the outermost field (bottom right). 
The star with ECa = 11 and M/ = —2 in the field situated at 5° has a well measured ECa 
and radial velocity = 264 km s~^ (the mean of the sample is 278 ± 20 km s~^), 
so it seems to be an LMC RGB star with unusually large metallicity. However, a reliable 
metallicity cannot be obtained with the present calibration, which is only valid for [Fe/H] < 
+0.47. The extrapolation of our calibration gives [Fe/H] = +0.72 ± 0.4 for this star. 

The metallicity distribution of each field is shown in Figure [5] and will be discussed in ex- 
tenso in Section [721 In Figure [6] the observed stars have been plotted in the color-magnitude 
diagram employing different symbols as a function of their metallicity. It is important to 
highlight that there is no correlation between the metallicity of the stars and their position 



in the color-magnitude diagram, as was also found in other galaxies (e.g. iPont et al.l 12004 



Cole et al.ll2005l : iKoch et al.ll2006l ). This shows once more that the metallicities derived from 



the position of the stars in the RGB are unreliable in systems with complex SFHs. 



6. Determination of stellar ages 



For a given age, the position of the RGB on the color-magnitude diagram depends 
mainly on metallicity in the sense that more metal-rich stars are redder than more metal- 
poor ones. At the same time, for a fixed metallicity, older stars are also redder than younger 
ones. If we take into account the logical assumption that younger stars are also more metal- 
rich, age may partly counteract the effect of metallicity. The combination of both effects 
is the well known RGB age-metallicity degeneracy, which sets a limit on the amount of 
information that can be retrieved from the RGB using only photometric data. However, if 
stellar metallicities are obtained independently from another source, s uch as spect r oscop y, we 
should in princip le be able to break the age-metallicity degeneracy. iPont et al.l (120041 ) and 
Cole et al.l (120051 ) have derived stellar ages from isochrones for stars whose metallicity had 
been previously obtained from low-resolution spectroscopy. In our case, instead of directly 
comparing the position of the star in the CMDs with isochrones, we have obtained a relation 
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for the age as a function of color, magnitude and metallicity from a synthetic CMD. This 
allows to easily obtain ages for the large number of stars in our sample. 

^The reference synthetic CMD has been computed using IAC-STAr| jAparicio fc GallartI 



20041 ) ■ Its main input parameters are the SFR and the chemical enrichment law, as a func- 



tion of time. Since we are interested in a general relation, we have chosen a constant SFR 
for the whole range of ages (13 > (Age/Gyr) > 0) and a chemical enrichment law such that 
a star of any age can have any metallicity in the range —2.5 < [Fe/H] < +0.5. Both the 
constant SFR and the chemical enrichment law ensure that there is no age or metallicity 
bias. We have coniputed two synthetic CMDs, one using the Ba STI stellar evolution library 
(jPietrinferni et al.ll2004l ) and another using the Padova library (iGirardi et al.ll2002l ) respec- 
tively. The ranges of ages and met allicities covered, as well as other i mportant features of 
each library are listed in Table 1 of iGallart. Zoccali fc Aparicio I ( 120051 ). 

Since the uncertainty of the age assigned by stellar evolution models to bright AGB 
stars is large, we will only compute ages for RGB and AGB stars fainter than the RGB tip 
(M/ ~-4 in Figure [T]) . We have chosen the synthetic stars in the box below the tip of the 
RGB shown in Figure [H We have computed a polynomial relationship which gives the age as 
a function of [Fe/H] = log(Z/Z0), (V-I) and My. In order to minimize the a and to improve 
the correlation coefficient of the relation, different linear, quadratic and cubic terms of each 
observed magnitude have been added. When the addition of a term did not improve the 
relation we rejected it. The best combination of these parameters is in the form of Equation 
[2] and the values obtained for each stellar evolution library are shown in Table HI Note that 
each stellar evolution library has its own age scale, and some differences might exist among 
them. 



log{agt) = a + 6(V - /) + cMy + d[Fe/H] + f{V - if + g[Fe/Hf + h{V - if (2) 

The uncertainties of each term on Equation [2] are given in Table HI However, the way 
in which they propagate into the error of the age is complex. For this reason a Monte Carlo 
test has been performed to estimate age errors. The synthetic stellar population computed 
from BaSTI stellar evolution models has been used to this purpose. That based on the 
Padua library would have produced similar results. The test consists in computing, for 
each synthetic star, several age values for stochastically varying [Fe/H], (V — I) and My 
according to a gaussian probability distribution of the corresponding a {<7[Fe/H] ~0.15 dex; 



•^available on the Web at 



http://iac-star.iac.es 
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(^{v-i) ~0.001 and ~ 0.001). The a value of the obtained ages provide an estimation 
of the age error when Equation [2] is used. The values obtained for different age intervals are 
shown in Figure [71 The error increases for older ages. 

We have also investigated the accuracy with which the age of a real stellar population is 
measured. The stars in open and globular clusters in Paper II, which ages are independently 
kno wn, are used for this purpose. W e adopt the ages estimated by ISalaris fc WeissI (120021 ) 
and lSalaris. Weiss, fc Percivall (120041 ). which are in a common scale, for a globular and open 
clusters, respectively. They are listed in Table 1 of Paper II. 

The age of each cluster star has been estimated from Equation [2l assuming the distance 
modulus and reddening listed in Paper II. The metallicity of each star has been calculated 
following the same procedure as for the LMC stars (see Section ED . The age of each cluster 
has been obtained as the mean age derived from all the stars in each cluster and the quoted 
uncertainty is the standard deviation of this mean. Ages derived for each cluster from the 
BaSTI and Padova relationships in Table HI versus the reference values for each cluster, have 
been plotted in Figure [HI The relationships saturate for ages older than 10 Gyr because 
differences of 1 or 2 Gyr produce only negligible variations of position in the CMD. In both 
cases, the age of clusters younger than 10 Gyr are well reproduced. Differences between 
the age derived with each model and the reference value are within the error bars, with the 
exception of cluster NGC 6819 for the BaSTI relationship. The large errorbar of the youngest 
cluster, NGC 6705, is due to the fact that the stars observed in this cluster are fainter that the 
synthetic stars used to derive the relationships. Therefore, we are extrapolating the relations 
to obtain the age of this cluster. As both models produce simila r values, we adopted fo r 
simplicity the relationships derived from the BaSTI stellar library (jPietrinferni et al.ll2004l ). 



7. Analysis 

7.1. Stellar kinematics of the LMC 

In Figured we h ave plotted the mean of the velocity distribution for each field (including 
the ICole et al.l [20051 . bar field) as a function of its Galactic longitude. The mean velocity 
changes from field to fi eld due to th e disk rotation of the LMC. The radial velocities of 
carbon stars derived by [Kunkel et al.[ (Il997l ) have also been plotted. Our data are in good 
agreement with their result, which is consistent with the presence of a rotational disk. 

If a classical Milky Way-like halo existed in the LMC (i.e. old, metal-poor and with high 
velocity dispersion), a dependency between the velocity dispersion and metallicity/age of the 
stars in our sample might be observed. The first evidence of the presence of a kinetically 
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hot spheroidal population was reported by lHughes. Wood, fc ReidI (Il99ll ). They found that 



the LMC long period variables, related to an old stella r population, have a high velocity 



dispersion (33 km s~^), with a low rotational component. iMinniti et al.l (120031 ) measured the 
kinematics of 43 RR Lyrae stars in the inner regions of the LMC and found that the velocity 
dispersion of these stars is 53 ± 10 km s~^, which they associated with the presence of a 
kinematically hot halo populated by old metal-poor stars. Some studies of the intermediate- 



Schommer et al. 


1992; 


Graff et al. 


2000). In fact 


Graff et al. 



disk population containing 20% of stars with a velocity dispersion of 8 km s ^, and an old 
disk population containing the remaining stars with a velocity dispersion of 22 km s~^. 

We can check the presence of a hot halo with our sample stars. Assuming that stars 
with similar metallicities are of similar ages (see below), a possible dependency of velocity 
dispersion with metallicity could indicate whether different stellar populations have different 
kinematics. 

The procedure has been the follow up. The mean velocity in each field has been sub- 
tracted from the radial velocity of each star to eliminate the rotational component of the 
disk. The total metallicity range has been divided into three bins. The velocity dispersion 
of stars in each metallicity bin is listed in Table [5] and plotted in Figure [TOl The velocity 
dis persion for old and meta l -poor stars, 26.4 km s~^, is slightly smaller than the value found 
by iHughes. Wood, fc ReidI (Il99ll ) for the old long period variables ((T=33 km s~^) and sig- 



nifica n tly smaller t han t he dispersion found in RR Lyrae stars {a ~53 km s~^. iMinniti et al. 
20031 ) . ICole et al.l (120051 ) found a velocity dispersion of a = 40.8 km s~^ for the most metal- 
poor stars in the bar, which is also higher than the value found here. The metal-poor stars 
in our sample seem to be members of the thick disk instead of the halo. In short, from our 
data we have found no clear evidence for the presence of a kinematically hot halo populated 
by old stars. 



7.2. Metallicity distribution 

The metallicity distribution for each field is shown in Figure |5l We have fitted a Gaussian 
to each one in order to obtain its mean value and dispersion (Table E]). The mean value is 
constant at [Fe/H]~-0.5 dex up to 8°, where the mean metallicity decreases by about a 
factor of two ([Fe/H]=-0.8 dex). All metallicity distributions are similar. They show a clear 
peak with a tail toward low metallicities. We complement our sample with the results by 
Cole et al. J2005h . who derived stellar metallicities in the bar in a similar way as here. As 
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demonstrated in Paper II, their CaT index is equivalent to ours. Using their measurements of 
the ECa we can calculate the metallicities from Equation [TJ We have also fitted a Gaussian 
to the resulting bar metallicity distribution to obtain its mean value. On average, the bar is 
slightly more metal-rich than the inne disk [Fe/H]~-0.4 dex). 



7.3. Age— metallicity relationships 

To understand the nature of the observed metallicity distributions, and to join insight 
on SFR and the chemical enrichment history of the LMC fields, we have estimated the age of 
each star in our sample following the procedure described in Section [HI These individual age 
determinations have a much larger uncertainty than the metallicity calculations. However, 
they are still useful because we are interested in the general trend rather than in obtaining 
precise values for individual stars. 

We assumed that the oldest stars h ave the same age as th e oldest cluster in our galaxy. 



for which we adop t 12.9 Gyr (NGC 64 26 jSalaris fc Weissll2002l ). This agrees with a 13.7±0.2 



Gyr old Universe (jSpergel et al. 2003 ) where the first stars formed about 1 Gyr after the Big 



Bang. Another important p oint is the youngest age t hat we can find in our sample. According 



to stellar evolution models (jPietrinferni et al.ll2004l ). we do not expect to find stars younger 



than 0.8 Gyr in the region of the RGB where we selected our objects. However, application 
of Equation [2] may result in ages younger than this value. To avoid this contradiction, and 
taking into account that the age determination uncertainty for the younger stars is about 1 
Gyr, we assigned an age of 0.8 Gyr to those stars for which Equation [2] gives younger values. 
Finally, as the relationships used to estimate the age have been calculated for stars below 
the tip of the RGB, only stars fainter than this point have been used. 

The age of each star versus its metallicity (i.e. the AMR) has been plotted in Figure 
[11] for each field. The many stars with the same age that appear both at the old (12.9 Gyr) 
and young (0.8 Gyr) limits are the consequence of the boundary conditions imposed on the 
ages, and the exact values should not be taken at face value. The age error in each age 
interval, computed in Section [6l is indicated in the top panel. The age distribution for each 
field has been plotted in inset panels. The histogram is the age distribution without taking 
into account the age uncertainty. The solid line is the same age distribution computed by 
taking into account the age error. To do this, a Gaussian probability distribution is used to 
represent the age of each star. The mean and a of the Gaussian are the age obtained for the 
star and its error, respectively. The area of each Gaussian is unity. In stars near the edges, 
the wings of the distributions may extend further than the age physical limits. In such cases, 
we truncated the wing and rescaled the rest of the distribution such that the area remains 
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unity. 

Stars brighter than M/=-3.5 were not used in Paper II to obtain the metallicity calibra- 
tion of SCa (Equation [T]). However, we observed stars brighter than M7=-3.5 in the LMC 
Therefore, extrapolation of Equation [1] has been used to obtain the metallicity of these stars, 
but only for those below the tip of the RGB. It is necessary to check whether this has intro- 
duced any bias in the obtained AMR. Paper II demonstrated that the sequences described 
by cluster stars used for the calibration in the Mj-UCa plane are not exactly linear and have 
a quadratic component. In the interval M7=0 - -3.5 the difference between the quadratic 
and linear behaviors are negligible. However, for brighter stars, the deviation from the linear 
behavior might be important, resulting in a possible underestimation of the metallicity of 
these stars. To check this point, in Figure [12] we have plotted as filled squares stars with 
magnitudes within the range of the calibration. Open circles are stars whose metallicities 
have been obtained by extrapolating the relationship (-3.5> Mj >-4). The mean metallic- 
ity and its dispersion for several age intervals have been also plotted for both cases. The 
differences between both groups are smaller than the observed dispersion and therefore we 
conclude that no strong bias is produced as a consequence of extrapolating Equation [1] for 
bright stars. 

Figure [TT] shows that the AMR is, within the uncertainties, very similar for all fields. As 
expected, the most metal-poor stars in each field are also the oldest ones. A rapid chemical 
enrichment at a very early epoch is followed by a period of very slow metallicity evolution 
until around 3 Gyr ago, when the galaxy started another period of chemical enrichment 
that is still ongoing. Furthermore, the age histograms for the three innermost fields are 
similar, although the total number of stars decreases when we move away from the centre. 
The outermost field has a lower fraction of "young" (1-4 Gyr) intermediate-age stars. This 
indicates that its lower mean metallicity is related to the lower fraction of intermediate- age, 
more metal-rich stars rather than to a different chemical enrichment history (for example, a 
slower metal enrichment). 

With the aim of obtaining a global AMR for the disk, in the following we will quan- 
titatively address the question of whether the AMRs of all our disk fields are statistically 
the same. We have calculated the mean metallicity of the stars and its dispersion in six 
age intervals, for each field and for the combination of the four. The values obtained are 
listed in Table [71 To compare the values obtained in each field with those for the combined 

sample, we performed a x test such that x = z2i=i 2^ — j where af is the sum 

of the uncertainties squared in the age bin i of the field and the combined AMR. The last 
column in Table [7] shows the values xl = X^/ \/5- From these values we may conclude that 
the AMR of each field is equivalent to that of the combined sample to a probability of more 
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than 99 per cent. 



Now we will compare the global AMR obtained for the disk with that observed in the 
bar. The bar stellar ages have been obtained using the relation calculated in Section O The 
age versus metallicity for each star in the bar (left) and disk (right) have been plotted in 
Figure [131 Notice that with a different method and different stellar evolution models, the 



AMR obtained here for the bar is quite similar to that obtained by I Cole et al.l (120051 ). We 
performed a fest as before between the galaxy disk and bar AMRs. They are the same to 
within a probability of 90 per cent. However, a visual comparison would indicate that, while 
in the disk the metallicity has risen monotonically over the last few Gyr, this is not the case 
for the bar. However, from the comparison of the mean metallicities and the errorbars for 
stars in this youngest bin (1-3 Gyr) it appears that this difference is statistically significant. 

The global disk AMR is qualitativel y similar to that of the LMC cluster system (e.g. 
Olszewski et al.l Il99ll : iDirsch et al.ll2000l ). except for the lack of intermediate-age clusters 
(Figure UM . Clusters show a rapid enrichment phase around 10 Gyr ago which is also 
observed in the disk and the bar. The second period of chemical enrichment in the last 3 
Gyr is also observed both in the field and in the clusters. The field intermedi ate-age stars may 



have contributed to the chemical enrichment at recent times. Note that the Grocholski et al. 



(120061 ) sample has a maximum metallicity [Fe/H]=-0.6, i.e. it does not seem to participate 
on the high metallicity tail at young ages that we and the other authors analysing clusters 
samples find. In the three works, there ar e a number of young clusters with low metallicity 
for their age, which iBekki fc Chibal (120071 ) associate with metal-poor gas stripped from the 
SMC due to tidal interaction between the SMC, LMC and the Galaxy over the last 2 Gyr. 
A few stars also have low metallicity for their age in our sample, but they represent a minor 
contribution to the total. 



7.4. The Star Formation History 

More information can be retrieved from the AMRs. In particular we can also derive an 
approximate SFH of the LMC if we can establish that the observed stars are representative 
of the total population. The number of stars in each age bin can be transformed into 
stellar total mass, after accounting for the number of stars with a given age that have 
alrea dy died. To evaluate this correction we have computed a synthetic CMD using lAC- 



star (lAparicio fc GallartI |2004| ) . We have assumed a constant SFR and used the relations 
derived for the chemical enrichment. As we have only observed a fraction of the total stars 
in the region of the RGB, we have rescaled the result to the total number of stars in this 
region. 
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In order to check whether the procedure used to select the stars observed spectroscop- 
ically introduces any bias in the computed SFH, we have calculated a synthetic CMD with 
a known chemical enrichment history and SFR. It was computed such that the number of 
stars in the selection region was the same as in the field at 3°. The input SFR as a function 
of age, ipit), is plotted as a solid line in Figure [151 To check that the selection criteria did 
not introduce biases, we have selected objects from this synthetic CMD in the same way as 
from the LMC fields (see Figured]). We assigned random positions to each synthetic star 
and, using the HYDRA configuration software, we computed 20 test configurations. The 
number of objects selected in each configuration, about 115, is similar to the observed stars 
in the LMC fields. For each test configuration, we calculated ip{t) following the procedure 
described above. The mean of the 20 total tests is the dashed line in Figure [151 The error 
bars are the standard deviation of the 20 tests. We repeated the same test, but obtaining 
the ages of each synthetic star using Equation [21 The result is the dotted line in Figure [13 
As we can see in this Figure, the recovered ip{t) are consistent with the real one. In half of 
the bins the differences are within the error bars and in the rest they are within 3a. 

The ilj{t) derived from our four fields are shown in Figure [T6l (soJid line), which has been 
computed from the age distribution, taking into account the age error (see Figure [7]). The 
uncertainty in the computed ip{t) has been estimated as the square root of the SFR value for 
each age {dotted lines). However, values of the SFR for different ages are not independent 
from each other because the integral for the full age interval is a boundary condition of the 
solution. In other words, fluctuations above the best solution should be compensated by 
others below it, and a solution close to, for example, the upper (or the lower) dashed line is 
very unlike. As a comparison, the ip{t) computed without taking into account the age error 
has also been plotted in Figure [161 AH flelds have a flrst episode of star formation more than 
10 Gyrs ago. Then, ilj{t) decreases until ~5 Gyr ago, when it rises again in the inner regions 
of the LMC (r<6°). This enhancement of ipit) is not observed in the outermost fleld. 

We followed the same procedure to calculate iplt) for the bar stars. The bar iplt) is 
shown in Figure [T71 In this case we do not know the total number of objects in the region 
where the observed stars were selected. For this reason, we have plotted the percentage, over 
the total, that the SFR represents in each age bin. In the right panel of the same Figure, the 
iplt) of the disk, obtained as a combination of the four flelds in our sample, is shown. The bar 



^(t) i s similar to that derived i n small flelds from HST deep photometry (IHoltzman et al. 



19991 : ISmecker-Hane et al.ll2002l ). as indeed is expected since the HST Helds are in the same 



region in which the RGB stars were observed. Both in the bar and in the disk, the galaxy 
started forming stars more than 10 Gyr ago, although in the bar the initial enhancement 
of ip{t) is not observed. For intermediate ages, the SFR was low until about 3 Gyr ago, 
when another increase in the SFR, which is particularly intense in the bar, is observed. This 
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is consistent with previous investigations (e.g. iHardy et al.lll984l : ISmecker-Hane et al.ll2002l ) 
and with the SFR derived from the clusters. The cluster age gap coincides with the age 
interval in which star formation has been less efficient. 



7.5. A Chemical Evolution Model for the LMC 



In what follows we use the AMR to obtain information on the physical parameters 
governing the chemical evolution of the LMC. As a first approximation, we will try to 
reproduce the derived relations with simple models. 



Under the assumption of insta n taneo us recycling approximation, following iTinsleyl (119801 ) 
and iPeimbert. Colin fc Sarmientd (119941 ). the heavy-element mass fraction in the ISM, Z, 
evolves via: 



^ y{l - R)ij + {Zf - Z)fj 
^dfi -(l-/?)^ + (/,-/o)(l-/i) 

where // and fo are the infiow and outfiow rates respectively, Zj is the metallicity of 
the infall gas, y denotes the yield, R is the mass fraction returned to the interstellar medium 
by a generation of stars, relative to the mass locked in stars and stellar remnants in that 
generation, ip is the SFR and /x is defined as /i = Mg/Mh, where M^ is the gas mass and M^ 
is the total baryonic mass, i.e., the mass participating in the chemical evolution process. In 
general, Z, /x, ip, fj and fo are explicit functions of time, while R and y may change according 
to characteristics of the stellar population like metallicity and initial mass function, but are 
usually assumed to be constant for a given population. 

When there are no gas fiows in and out of the system (// = /o = 0), the model is 
called a closed-box model and the solution to Equation [3] provides the variation of Z{t) as a 
function of fj,(t) via 



Z{t) = Zi + y\iafi{t) 



-1 



(4) 



In an infall scenario, in which gas fiows into the system at a rate proportional to the 
amount of star formation, fj = a{l — R)ip, a being a free parameter, we can integrate 
Equation [3] for a ^ 1, assuming that fo = 0, and obtaining: 



Z{t) = Z, 



y + Zfa 



a 



a 



a 



fi{t) 



-a/{l-a)' 



(5) 



-17- 



For a=l, the solution of Equation [3] is 



Zi + y 



(6) 



In an outflow scenario, in which gas scapes from the system at a rate fo = A(l — R)ip, 
A being a free parameter, we obtain from Equation [3] assuming // = 0, and for (A 7^ 1): 



Z{t) = Zi + 



y 



A 



In 



A 



A 



(7) 



Finally, we will consider a model combination of inflow and outflow, and assume that 
the system gas flow is proportional to the mass that has taken part in the star formation: 
/ = (a — A)(l — R)ip, where we can deflne (3 = a — X. In this case, we can integrate Equation 
[3] for P 7^1, to obtain 



Z{t) = Z,+ 



y + Zfa 



a 



/5 



1-/3 



-/3/{l-/3)' 



(8) 



In summary, the chemical evolution of the galaxy is given by the yield ?/, the initial 
metallicity of the gas Zj, the fraction of gas mass /i(t), which is a function of ip{t), and i?, 
the mass fraction returned to the interstellar medium. In the case of inflow and/or outflow, 
also by the additional parameters a and A. 

We can now apply the former relations to the LMC In the previous section we have 
derived il){t) for the LMC bar and disk, the latter being obtained from the combination of the 
four flelds in our sample, /i is computed explicitly as a function of time from the balance 
of gas flowing to and from the system, and using as a boundary condition the cu rrent LMC 



gas fr action. In the case of the LMC disk, we have assumed the value derived by iKim et al. 



(119981 ) from observations of H I. These authors derived Mg = 5.2 x lO^M© and estimated 



that the total mass of the disk is Mb = 2.5 x IO^Mq. With these values, we flnd the current 
gas fract ion in the disk, = 0. 21. For the bar we have assumed the value of n = 0.08 



Pagel fc Tautvaisiend ( 



given by IWesterlund et al.l (119901) . For the yield, we have assumed the value obtain ed by 



I995I ) for the s o lar n eighborhood, y = 0.014. Assuming the IScalo 



(1l986l ) initial mass function, iMaederl (1l993l ) calculated a returned fraction associated to 
this yield of i? = 0.44. Putting all these ingredients together, we will now explore several 
scenarios in order to reproduce the AMR of the disk and the bar. 



In Figure [18] we have plotted the age and metallicity of each star in the disk, together 
with the mean metallicity and its dispersion, for stars in each age bin. In this flgure, chemical 
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evolution models have been superposed for different scenarios and parameter assumptions 
(see the figure caption). In the case of the infall models {green dashed lines), we have 
assumed a zero-metallicity infalling gas, in order to obtain the current metallicity. They do 
not reproduce the general behavior of the observed AMR. The disk AMR is well reproduced 
by outflow models {blue dot-dashed lines) with a relatively large range of A. Finally, the 
combination of inflow and outflow {pink and cyan long-short dashed lines) models also 
reproduce the metallicity tendency (in this case the model that best matches the observed 
data has a = 0.2 and A = 0.05). The previous models have been computed assuming the 
yield observed in the solar neighborhood which is probably too large for the LMC With a 
smaller yield, y = 0.008 {red thick solid line), the observed AMR can be re produced by a 



closed-box model. This yield would correspond to stars with Z = 0.001 (see lMaederlll993 



for details), which is slightly low compared with the mea n LMC metallicity. Bursting and 



smooth models computed by iPagel fc Tautvaisiend (Il998l ) have been plotted as comparison 



{brown lines). The bursting model was computed assuming two burst of star formation, one 
about 12 Gyr ago and a strong one about 3 Gyr ago, and some fraction of infall and outflow. 
Note that their predicted values agree with our mean metallicity for each age bin, within 
the errorbars, specially at old and intermediate-age. It also qualitatively agrees with the 
episode of faster chemical evolution with started around ~2 Gyr ago. The smooth model was 
computed assuming a constant star formation rate and does not match with the observations. 

In the case of the bar (Figure [T9|) . the very slowly rising metallicity, almost constant over 
the last few Gyr, is best reproduced by the combination of models with inflow and outflow 
with infall parameter a = 1.2 and A = 0.4-0.6 {pink long-short dashed lines). Infall models 
{green dashed lines) with a = 0.4-0.6 marginally reproduce the observed trend, but they 
predict higher metallicities than observed in the last few Gyr, unless a smaller yield would 
be assumed. Finally, outflow and closed -box models predict a rise o f metallicity in the last 



few Gyr more steep than observed. The iPagel fc Tautvaisiend (119981 ) bursting model agrees 



worse with the bar AMR than with that of the disk, specially in the youngest age bins. 



8. Conclusions 

Using infrared spectra in the CaT region, we have obtained metallicities and radial 
velocities for a sample of stars in four LMC fields. Metallicities have been calculated using 
the relationships between the equivalent width of the CaT lines, SCa, and metallicity derived 
in Paper II. In addition, we have estimated the age of each star using a relationship derived 
from a synthetic CMD which, from the color, magnitude and metallicity of a star, allow us 
to estimate its age. The main results of this paper are: 
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• The velocity distribution observed in each field agrees with the rotational thick disk 
kinematics of the LMC. The velocity dispersion is slightly larger for the most metal- 
poor stars. However, the values obtained do not indicate the presence of a kinematically 
hot halo. 

• The metallicity distribution of each field has a well defined peak with a tail toward low 
metaUicities. The mean metaUicity is constant until the field at 6° ([Fe/H]~-0.5 dex), 
and is a factor two more metal-poor for the outermost field ([Fe/H]~-0.8 dex). 

• The AMR observed in each disk field is compatible with a single global relationship for 
the disk. We conclude that the outermost field is more metal-poor on average because 
it contains a lower fraction of relatively young stars (age<5 Gyr), which are also more 
metal-rich. 

• The disk AMR shows a prompt initial chemical enrichment. Subsequently, the metal- 
licity increased very slowly until about 3 Gyr ago, when the rate of metal enrichment 
increased again. This AMR is similar to that of the cluster system, except for the lack 
of clusters with ages between 3 and 10 Gyr. The recent fast enrichment observed in 
the disk and in the cluster system is not observed in the bar. 

• The ip{t) of the three innermost fields show a first episode of star formation until about 
10 Gyr ago, followed by a period with a low SFR until ~5 Gyr ago, when the SFR 
increases, and reaches its highest values ~2-3 Gyr ago. The outermost field does not 
show the recent increase of SFR. The second main episode is also observed, and is 
more prominent in the bar, where an increased SFR at old ages is not observed. The 
lower SFR between 5 and 10 Gyr ago is probably related to the age-gap observed in 
the clusters. 

• Under the assumption of a solar yield, the disk AMR is well reproduced either by a 
chemical evolution model with outfiow with A between 1 and 2, so the disk losses the 
same amount of gas that has taken part in the star formation, or by models combining 
infall and outfiow with q;=0.2 and A=0.05, which means that the galaxy is almost a 
closed-box system. With a smaller yield, the AMR could also be reproduced with a 
closed-box model. The bar AMR is well reproduced by models with a combination 
of infiow and outflow with q;=1.2 and A between 0.4 and 0.5. This suggests that 
the amount of infalling gas was larger than the amount that participated in the star 
formation in the bar, and also that the amount of gas that escaped the bar was 50% 
of the total that participated in star formation. 
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9. Discussion 



The main result in tliis paper is that all our fields, covering galactocentric radius from 
3° to 8° from the center (2.7 to 7.2 kpc) share a very similar AMR. The mean metallicity is 
very similar in the three innermost fields ([Fe/H] ~-0.5 dex), and it is a factor of ~ 2 smaller 
in the outermost field ([Fe/H] ~-0.8 dex). Because the AMRs are the same, this has to be 
related with the lower fraction of relatively young stars (age < 5 Gyr), which are also more 
metal-rich, in the outer part. We find, therefore, a change in the age composition of the disk 
population beyond a certain radius (~ 6 kpc; note however that in this study we are only 
sensitive to populations older than ~ 0.8 Gyr, and that differences among fields for younger 
ages could also exist), while the chemical enrichment history seems to be shared by all fields. 

The comparison of the AMR observed in the disk with simple chemical evolution models 
suggests that the LMC is mo st likely losing some of its gas. This is in agreement with the 
work by iNidever et al.l (120071 ) which suggests that the main contribution to the Magellanic 
Stream comes from the LMC instead of from the SMC as it was believed until now (e.g. 



Putman et al.ll2003l ). Alternatively, the LMC could be an almost closed-box system with 
small gas exch a nges (a =0.2, A=0.05). This could be in agreement with the models by 
Bekki fc Chibal (120071 ) which suggest that the LMC has received metal-poor gas from the 
SMC in the last 2 Gyr. In our field we observed some stars with a low metallicity to which 
we assign an age younger than 2 Gyr. However, the uncertainty in the age determination is 
large. 



We have also obtained the AMR of the LMC bar from the data by ICole et al.l (120051 ). 
The bar AMR differs from the one of the disk in the last 5 Gyr: while in the disk the 
metallicity has increased in this time, in the bar it has remained approximately constant. 
This feature is best reproduced by models combination of outflow and a relatively large infall 
of pristine gas. Models with a smaller yield, but with an infall of previously enriched gas also 
reproduce the observed bar AMR. This would be in agreement wit h the prediction that a 
typic al bar instability pushes the gas of the disk towards the center ( ISellwood fc Wilkinson 
19931 ). so the infall gas is expected to be previously pre-enrichmen t by the disk c h emica l 
evolution. This would be also the case in the scenario suggested by lBekki &: Chibal (120051 ). 
in which the bar would have formed from disk material as a consequence of tidal interactions 
between the LMC, the SMC and the Milky Way about 5 Gyr ago, the moment when the 
bar AMR differs from that of the disk. Also, a bar is expected to destroy any metallicity 
gradient within a certain radius, as it is observed. 



Gallart et al.l (l2004j ) derived the LMC surface brightness proflle using deep resolved star 



photometry of the four flelds in the current spectroscopic study, and found that it remains 
exponential to a radius of 8° (^ 7 kpc), with no evidence of disk truncation. Combining this 
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information with that on the deep CMD of the outermost field, which contains a large frac- 
tion of intermediate-age stars, they concluded that the LMC disk extends (and dominates 
over a possible halo) at a distance of at least 7 kpc from its center. In the present study, the 
kinematics of the stars in these four fields confirm this conclusion: the velocity dispersion 
of all four fields is similar, around 20 km/sec. If the stars are binned by metallicity, the 
velocity dispersion of the most metal-poor bin is slightly larger than that of the metal-richer 
bins (cry ~ 25 km/sec as opposed to ~ 20 km/sec) but still not large enough to indicate 
the presence of a halo, even one formed as a consequence of the interaction with the Galaxy. 
In this case, iBekki et al.l (120041 ) predicted a velocity dispersion ~ 40 km/s at a distance 
of 7.5 kpc from the LMC center (similar to the dista nce of our outermost field). Other au- 



thors have also failed to find a kinematically hot halo (IFreeman. lUingworth. fc Oemlerlll983 



Schommer et al.l Il992l : I Graff et al.l l2000l : IZhao et al.l |2003| ) . The fir st evidence of the pres- 



ence o f a kinetically hot, old spheroidal population was reported by lHughes. Wood, fc Reid 
(Il99ll ) us ing a sample o f long period variables, which are related to an old stellar population. 
Recently, iMinniti et al.l (120031 ) observed spectroscopically a sample of 43 bar RR Lyrae stars, 
and obtained a large velocity dispersion (53 ± 10 km s~^) for them. It is possible that our 
failure (and that of other authors) to find evidence of a hot stellar halo is related with a low 
contrast of the halo population with respect to the disk one, even at large galactocentric 
radius as our outermost field. 
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Fig. 1. — CMD of the RGB region of the field situated at 6° from the LMC center. The 
regions used to select stars to be observed spectroscopically are shown. 
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Fig. 2. — Example of sky subtraction. Top: raw spectrum. Middle: sky spectrum obtained 
from a combination of fibers located on the sky. Bottom: The final spectrum after sky 
subtraction. 
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Fig. 3. — Radial velocity distributions of our four LMC fields. The mean value of each 
distribution is shown in each panel, and its position is marked with an arrow. Note that 
there are two peaks in the field situated at 3°. 
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Fig. 4. — Position in the Mj-SCa plane of LMC stars for which membership has been 
confirmed from their radial velocity. The typical SCa error is shown in the bottom right 
corner of each panel. Isometallicity lines have been plotted for reference. The solid part of 
each line is the magnitude interval covered by the cluster stars used for the calibration (see 
Paper II). The dashed part is the region in which the calibration is extrapolated. 
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Fig. 5. — Metallicity distributions for the four fields in our sample. A Gaussian has been 
fitted to each distribution in order to obtain its mean and dispersion. The values obtained 
are shown in each panel. 
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Fig. 6. — Position of the observed star in the color- magnitude diagram. Different metaUicity 
intervals are represented by different symbols, as indicated in the label. Solid lines are the 
fiducial sequences of four globular clusters (from right to left: NGC 104: [Fe/H] = —0.78; 
NGC 5904: [Fe/H] =-1.12; NGC 5272: [Fe/H] = -1.34 and NGC 6341: [Fe/H] = -2.16). 
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Fig. 7. — Uncertainty in the age calculated using Equation [2] obtained from the BaSTI stellar 
evolution models. See text for details. 
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Fig. 8. — Ages derived for the cluster sample in Paper II, using Equation [2] and BaSTI and 
Padova models versus the reference ages. The solid lines are the one-to-one relation. 
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Fig. 9. — Radial velocities for carbon stars from iKunkel et al.l (119971 ) in the band —36° < 
b < —33°. The triangles are the ca r bon s tars in our fields at 5° and 8°. The filled square 
represents the data from lCole et al.l (120051 ) for the bar. The filled dots are our data listed in 
Table O The error bars represent the dispersions of the velocity distributions. 
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Fig. 10. — Radial velocity dispersion of the disk stars as a function of metallicity. Although 
the velocity dispersion increases from metal-rich to metal-poor stars, this increment is not 
clear enough to conclude that a metal-poor stellar halo is present. 
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Fig. 11. — Age-metallicity relationships for the four LMC fields in our sample. Inset panels 
show the age distribution computed taking into account {solid line) and not {histogram) the 
age determination uncertainties. The top panel show the age error in each interval. The left 
panels show the metallicity error. 
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Fig. 12. — AMR of the field situated at 5° for different star groups. Squares: objects 
witfi magnitudes in the range of the M/-SCa cahbration obtained in Paper II; open circles: 
objects in the same field with magnitudes brighter than the cluster calibration stars. The 
corresponding average metallicity and dispersion in different age bins are also plotted. 
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Fig. 13. — Same as Figure [TT] but for the bar (left) and the disk fields combined {right). 
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Fig. 14.— L MC clusters AMR bv loiszewski et all (|l99ll . squares); IPirsch et all kom 



triangles) and iGrocholski et al.l (l2006l . hexagons). The mean metallicity in six age bins of 
our global disk has been plotted {filled points) . Note that the metallicity scales of each work 
may not be exactly the same. 
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Fig. 15. — Results of the test on the rehabihty of the SFH computed using our spectroscopic 
RGB sample. Solid line: injected SFH; dashed line: mean and sigma of the recovered SFH 
after the Hydra filter assignment simulation; dotted line: the same as dashed one but with 
the age obtained from Equation [21 See text for detail. 
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Fig. 16. — SFR, as a function of time, for our four fields. Note that the y-axis scale is different 
in the top and bottom panels. Histogram is the tp{t) computed from the age distribution 
without taking into account the age uncertainty. Solid line is the same but computed by 
taken into account the age error, while dotted lines are its uncertainty. 
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Fig. 17. — SFR as a function of time for the bar (left) and disk (right) derived from the 
spectroscopic sample. 
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Fig. 18. — Chemical evolution models overplotted on the AMR for the disk. Red solid lines 
are the closed-box models with y = 0.008, {thick line) and y = 0.014 {thin hne), respectively. 
Blue dot-dashed lines are outflow models for A = 0.6, 1.3, 2.0 and 2.5. Green dashed lines 
are infall models with a = 1, 0.9, 0.5 and 0.3. Pink long-short dashed lines are models 
with infall and outflow for a = 0.2, and A = 0.1 {upper) and 0.05 {lower). Finally, cyan 
short-long dashed line is a model with inflow and outflow for a = 0.3 and A = 0.2. All 
infall and outflow models assume y = 0.01 4. The brown solid and dashe d lines represent the 
bursting and smooth models computed by lPagel fc Tautvaisiene! (Il998l ). 
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Fig. 19. — As Figure [18] for the bar. Red solid line is a closed-box model. Blue dot-dashed 
lines are outflow models for A = 0.8 and 0.3. Green short-dashed lines are infall models with 
a = 1, 0.6 and 0.4. Pink long-short dashed lines are models with inflow and outflow for a = 
1.2 and A = 0.4, 0.5 and 0.6. All model s assume y=0.014. As in the previous flgure, brown 
lines represent the models computed by lPagel &: Tautvaisiene! (1l998l ). 
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Table 1. LMC RGB stars observed 



02000 <52ooo S Ca osca V I Vr{kms ^) cFv^{kms -"-) Comments 



05 


08 


53.92 


-66:49:55.1 


7.8 


0.8 


17.67 


16.36 


287.3 


0.7 




05 


08 


55.21 


-66:47:34.7 


8.7 


0.2 


16.66 


14.70 


240.9 


0.6 




05 


08 


53.89 


-67:02:50.0 


2.4 


0.2 


17.19 


14.25 


23.0 


0.6 


No member 


05 


08 


55.40 


-66:53:45.4 


8.9 


0.4 


16.79 


15.09 


292.3 


0.5 




05 


08 


58.58 


-66:45:49.2 


9.4 


0.2 


16.00 


13.99 


291.9 


0.5 





Note. — Table[T]is published in its entirety in the electronic edition of The Astronomical Journal. A portion 
is shown here for guidance regarding its form and content 



Table 2. Radial velocities and velocity dispersion for each field. 



Field 


<K)(fcms-i) 


(Jv{kms ^) 


Bar 


260 


24 


3° 


269 


27 


5° 


278 


20 


6° 


293 


15 


8° 


282 


14 



Table 3. Line and continuum bandpasses 



Line Bandpasses (A) 


Continuum bandpasses (A) 


8484-8513 


8474-8484 


8522-8562 


8563-8577 


8642-8682 


8619-8642 




8799-8725 




8776-8792 



Table 4. Coefficients of Equation [2l 



Model 


a 


b 


c 


d 


f 


g 


h 


f7 


BaSTI 
Padova 


2.57±0.10 
1.69±0.07 


9.72±0.15 
10.6±0.1 


0.70±0.003 
0.75±0.002 


-1.51±0.007 
-1.87±0.007 


-3.86±0.08 
-4.14±0.06 


-0.19±0.007 
-0.24±0.006 


0.49±0.01 
0.51±0.009 


0.38 
0.37 
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Table 5. Velocity dispersion for each metallicity bin. 



[Fe/H] 


Age (Gyr) 


N« 


{Vr){kms-^) 


(Ty {kms ^ ) 


>-0.5 


<3 


226 


0.1 


20.5 


-0.5 to -1 


2-10 


178 


1.9 


20.5 


<-l 


>10 


100 


1.8 


26.4 



Table 6. Mean values of metallicity distributions. 



Field 


([Fe/H]) 


^[Fe/H] 


Bar 


-0.39 


0.19 


3° 


-0.47 


0.31 


5° 


-0.50 


0.37 


6° 


-0.45 


0.31 


8° 


-0.79 


0.44 



Table 7. Average metallicity in six age bins. 



Field 


aFe/H]<i.5> 






([J='e/if]5.5_8.5> 


([Fe/if]8.5-ll> 


([Fe/H]>ii) 


xl 


Bar 


-0.27±0.13 


-0.32±0.16 


-0.63±0.26 


-0.66±0.30 


-1.02±0.41 


-1.25±0.56 




Disk* 


-0.17±0.21 


-0.39±0.15 


-0.60±0.16 


-0.71±0.18 


-0.80±0.27 


-1.25±0.41 




3° 


-0.16±0.17 


-0.40±0.14 


-0.58±0.15 


-0.67±0.20 


-0.69±0.07 


-1.13±0.43 


0.05 


5° 


-0.1(i±0.28 


-0.39±0.15 


-().63±0.18 


-().72±().16 


-().83±0.35 


-1.39±0.42 


0.02 


6° 


-0.15±0.16 


-0.38±0.16 


-0.60±0.16 


-0.75±0.20 


-0.83±0.14 


-1.34±0.24 


0.02 


8° 


-0.16±0.21 


-0.45±0.13 


-0.62±0.20 


-0.89±0.16 


-0.92±0.18 


-1.25±0.32 


0.16 



^Combination of the four fields. 



